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We demonstrate the use of a phase- only spatial light modulator for the measurement of transverse spatial 
distributions of coincidence counts between twin photon beams, in a fully automated fashion. This is 
accomplished by means of the polarization dependence of the modulator, which allows the conversion of a 
phase pattern into an amplitude pattern. We also present a correction procedure, that accounts for 
unwanted coincidence counts due to polarization decoherence effects. 



Spatial variables of photon pairs have proven to be a useful tool to investigate quantum entanglement and 
quantum information\ particularly due to their high-dimensional nature^"'*, capacity for quantum state 
engineering^"^, and application to quantum communication^'^ and quantum imaging^^"^^. Spatial correla- 
tions can be measured by digitizing the detection planes and obtaining a raster graphic of the correlations. This 
kind of measurement can be accomplished by raster- scanning single-pixel detectors in the transverse plane^'*"^^, 
using single-photon sensitive CCD cameras^^"^^, with other methods such as compressed sensing^^, or using 
interferometric techniques^^'^^, for instance. These methods involve reconstructing the marginal probability 
distributions associated to the spatial variables. 

The spatial light modulator (SLM) has revolutionized the types of operations that can be implemented on the 
spatial variables of the field. In general, the SLM applies a user-defined position- dependent phase S(x, y) on the 
incident wavefront. In this way, it is possible to engineer high-dimensional entangled states^^ and also implement 
a number of quantum logic gates^^ in the spatial variables. Particular interest has been payed to the orbital angular 
momentum degrees of freedom, where SLM's have been used to test entanglement^^ and Bell's inequality^^ as well 
as mutually unbiased bases in six dimensions'^, and to simulate stronger-than-quantum correlations'\ An 
interesting aspect of certain types of SLMs is the fact that they can be polarization-dependent, imprinting a phase 
on the field for-say-horizontal (H) polarization and doing nothing when the polarization is vertical (V). 
Considering both the polarization and spatial degrees of freedom, the action of the SLM can be described by 
the operator 



S = \V){y\®\+\H){H\®e' 



'S(x,y) 



(1) 



where I is the identity. Using this coupling, quantum logic gates between polarization and spatial parity qubits 
have been implemented'^, a qubit coupled to a chaotic quantum harmonic oscillator has been studied'', a single 
photon and classical channels have been multiplexed in a few-mode fiber'^, spatial moments have been measured 
directly'^, and an optical integration algorithm has been performed'^. In Ref. [37], quantum tomography was 
performed on the polarization qubit for constant homogeneous phases of the SLM, showing that the SLM 
implements S with a fidelity of about 92%. The imperfect fidelity is due primarily to a slight decoherence of 
the polarization induced by the SLM. This decoherence can have a detrimental effect on quantum logic gates, and 
needs to be corrected in some situations'^. 

Here we show that the coupling between the polarization and the spatial degrees of freedom (DOF) can be 
used to perform spatial raster-scanning of a wavefront with the SLM. In this procedure, a phase slit is 
scanned across the SLM. This allows one to characterize the spatial profile of the field on the SLM in a fully 
automated fashion. In our measurements we observe spatially inhomogenous background counts, resulting 
from imperfections in the SLM. This reduces the visibility or constrast of the reconstructed pattern. We then 
develop a model for the background counts that originate from the polarization decoherence of the SLM. All 
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Figure 1 | SLM phase calibration curves, a) Phase calibration curve for photon s. b) Phase calibration curve for photon /. 



parameters from the model can be experimentally determined. 
Using this model, we can account for the decoherence in such a 
fashion that the background counts are discarded, thus increasing 
the contrast of the coincidence distributions. 

Raster Scanning with the SLM 

Many experiments require the reconstruction of the two -photon 
coincidence distribution P{Xs, x^) as a function of the positions Xs 
and Xi of the signal (s) and idler (/) photons. Traditionally, this has 
been done by scanning detectors equipped with slit apertures in the 
transverse plane, using micrometers or stepper motors, as in Refs. 
[15,17]. More recently, digital micromirrors^ or single-photon 
sensitive CCD cameras^^"^° have been used. As an alternative 
method, we propose to exploit the polarization properties of the 
SLM. The method consists of producing a series of phase- slit 
images on the SLM, which are converted into amplitude slits 
through post- selection of the polarization state of light. This allows 
for fully automated scanning. To produce a phase pattern with the 
SLM, an 8 bit (0-255) greyscale image g{x, y) is programmed onto 
its LCD screen using a computer. An incoming field receives a 
phase y) = 2ng(x, y)l256. To perform a linear scan in the 
X-axis of the signal beam, we program the SLM with a sequence of 
N greyscale images J„ composed of phase slits {n = 0.. N — 1) of 
width A where the phase pattern is given by 



(n+l)A (m+l)A 



an{xs) = 



0 for nA<Xs<{n+l)A, 
n elsewhere. 



(2) 



In other words, we program a zero-phase slit of width A on the 
SLM in the region defined above, while all other regions have 
phase 71. A similar phase is applied to the idler beam /, so that 
we have bmixi), defined in the same way. We consider the two- 
photon state obtained from spontaneous parametric down-conver- 
sion (SPDC) in the paraxial and monochromatic approximations^: 



= J J dxsdxi\l/{xs,Xi)\xs. 



(3) 



where \xp +) refers to a single photon state in the position repre- 
sentation and with the diagonal polarization state: | + ) = ^ly/ V2^ 

+ Applying the SLM operation S, on both signal and 
idler beams and projecting onto the |+)5|+)/ joint polarization 
state, the coincidence detection probability is 



P{n,m) = J dxs j dxi\\jj{xs,Xi) 



(4) 



In this way, the entire N X N two -photon coincidence distribution 
can be reconstructed. The spatial resolution of this method 
depends upon the slit width A. On the one hand, for large slit 
width, the reconstructed distribution P(/i, m) is strongly discre- 
tized, which can cause problems when one wishes to estimate 
entanglement^^'^^. On the other hand, when the slit width A is 
narrow, unwanted background counts that appear due to imper- 
fections in the operation of the SLM become increasingly more 
relevant, as we will see explicitly below. In the next section we 
describe a model for the origin of these background counts and 
present a post-processing method to reduce them, greatly increas- 
ing the fidelity with the true coincidence distribution. 

Improved Scanning Method 

Programming a constant phase (j) on the SLM, we can perform a 

polarization rotation of the form | + + exp(/(/>)|H))y/\/2, 

on the incident field. Ideally, a subsequent polarization measurement 
using a polarizing beam splitter would project this state onto the two 
orthogonal polarization states |+) and |— ) with probabilities P+ = 
cos^((/>/2) and P_ = sin^((/>/2). Therefore, measurements of P+ and 
P- give us an interference curve whose phase is defined by the SLM. 
In Figs, la) and b) we show this kind of polarization interference as a 
function of the uniform greyscale value applied to the SLM. In these 
measurements, one photon was detected in the | + ) polarization state 
and the other in the | V) polarization state that is unsensitive to the 
SLM, so as not to interfere with the polarization oscillations of the 
other photon. In this way, these curves correspond to single-photon 
interference curves, and provide information about the noise counts 
that occur due to polarization incoherence. 
The visibility is defined as 



V-- 



I max ~l~ -^rrt 



(5) 



where I is the coincidence count rate. Ideally, I^ax ^ Imin = 0 and 
V= 1. However, due to imperfections of the SLM, the visibility is 
reduced. Fitting the curves with a sinusoidal function, we find 
= 0.947 ±0.004 and V/ = 0.87 ±0.01 for the signal and idler 
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photons, respectively. In principle, the visibilities should be equal. 
The difference between them is probably due to tiny phase fluctua- 
tions in each single pixel^°. This means that the average noise 
depends on the effective number of pixels interacting with the beam, 
and on the transverse field distribution. Therefore, differences in the 
diameter and intensity distribution of the two beams, may lead to 
different visibilities. 

Let us define p as the probability that the modulation is succesfully 
applied to a photon field. This means that, at the polarizing beam 
splitter, the photon is transmitted when the SLM is set to apply a 
phase zero and reflected when the SLM is set to apply a phase n. Then 
we have I^ax — pl ^min — (1 ~ p)^- When p = 1, we recover the 
ideal case. Substituting the expressions for 1^^^ and J^/„ into Eq. 5 
and rearranging, we find p={y +1) /2. 

Let us see how these interference curves provide information 
about the scanning measurements. The polarization measurement, 
together with the ideal phase slit, given by the function a^ixs, ys), 
defines an aperture described by a rectangle function with transmis- 
sion 1 in the region nA<x< {n + 1)A and 0 otherwise. In practice, 
due to the less than unity visibility, the effective aperture is a rectangle 
function with transmission/) in the region nA < x < (n + 1) A and (1 
— p) otherwise. As p can be determined from the visibility V, we can 
experimentally determine the value of p, and take the decoherence 
into account. 

Let us first describe our model for the simplified case where the 
SLM surface is divided in 3 slit regions for each light beam, as shown 
in Fig. 2. We denote Cy as the ideal (V= 1) number of coincidence 
counts when slit 1 is at position / and slit 2 is at position;. Here 1 and 
2 refer to signal and idler, respectively. We define ATy as the measured 
number of coincidence counts between these two regions. Ideally, for 
a perfect SLM and perfect polarization optics, we would detect coin- 
cidence counts only between photons incident on the phase slits, so 
that Nij = Cfj. However, due to the imperfections discussed above, we 
detect coincidence counts that originate from other regions of the 
SLM. The probability that a photon is detected when coming from a 
zero phase region, implementing the equivalent of a slit, is given by 
thep obtained from the visibility measurements, as explained above. 
Likewise (1 — p) is the probability that a photon is detected when 
coming from a region modulated with phase n. Then, we see a 
decrease in signal, since instead of detecting all ideal counts Cn we 
detect />5p/Cii < Cn of them. In addition, we detect unwanted counts 
coming from different combinations of Ti-phase regions of the SLM. 
For example, it is possible that photon s reflected from region 1 and 
photon / from region 2 or 3, contributing to a background term. 
Then, we expect a contribution from terms C12 and C13 that is pro- 
portional tops(l — pi). Similarly, terms like C22 and C23 appear with a 
proportionality constant (1 — ps){l — pi). Summing up all these 
events, we can relate the measured count rate Nn to the ideal count 
rates as 



iVii = PspiCn +ps{l -pi){Cu + C13) + (1 -ps)pi{C2i + C31) 

+ (1-P,)(1-P0(C22 + C23 + C32 + C33). 



(6) 
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Let us arrange the measured and ideal coincidence counts as column 
vectors: 



N = 



/Nn\ 
Nu 
Nn 
N21 
N22 

VN33/ 



and C = 



/c„\ 

C12 

Cl3 

C21 

C22 
C23 

C31 

C32 

VC33/ 



(7) 



We can write the measured coincidence counts N as a function of the 
ideal coincidence counts C as 



N = EC, 



(8) 



where the symmetric matrix 
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Finding the inverse of this matrix E"\ we can find the ideal coincid- 
ence counts as a function of the measured counts 



C = E-'N. 



(9) 



In the more general case of d slit positions, we follow the same 
procedure as outlined above, to obtain 

Nij = OCQj + P^Cmn+y'^Qn + 3^ Cmj, 



(10) 



and the matrices E and E ^ are described accordingly. 
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Figure 2 | Phase slits on the SLM for measurement of coincidences Nu 



Figure 3 | Experimental setup. Signal (blue) and idler (red) beams are 
directed onto the SLM. As can be seen in the image, each of the twin 
photons is incident on one side of the SLM. Acronyms are defined in the 
text. 
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Figure 4 | Experimental results when the detectors are scanned in the transverse plane, a) The marginal coincidence distribution Cs{x^. b) The marginal 
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Experiment 

Figure 3 shows a sketch of the experimental setup. Horizontally 
polarized frequency- degenerate twin photons at 650 nm are gen- 
erated by pumping a non-linear BBO (j5-Barium Borate) crystal 
with a He-Cd laser operating at >^ = 325 nm. The photons are 
directed onto a single SLM using mirrors and lenses. The SLM is 
the Pluto reflective phase-modulation-only model fabricated by 
Holoeye, with resolution 1920 X 1080 pixels and 8.0- fim pixel 
pitch. A half-wave plate is used right before the SLM to prepare 
the polarization state of the photons to | +). At the detection stage, 
each beam goes through another wave plate and a polarizing beam 
splitter to select the photons that reflected off of the intended 
region of the SLM. The detectors are single-photon avalanche 
photodiodes (APD) with large detection aperture. Lenses (focal 
lengths 70 cm and 15 cm) are used to map the tranverse position 
of the down -converted beams onto the SLM. For the measure- 
ments, horizontal phase slits are scanned vertically on each half of 
the SLM, corresponding to each beam. The total scan region for 
each photon is thus 1080 pixels. 



Results 

In order to characterize the spatial distribution of the coincidence 
counts in the absence of the SLM, we first performed the traditional 
procedure of scanning the detectors in the transverse plane and 
registering coincidence counts as a function of the detector position 
X. In these measurements the scanning detectors were equipped with 
a 20 -^m slit aperture and scanned in the vertical direction. The 
polarization of the photons was set to the vertical direction, so that 
the SLM had no effect on the spatial distribution. Figure 4a) shows 
results when detector is scanned and detector D/ is completely 
open (area- integrating). Fig. 4b) shows results when the detectors are 
scanned in opposite directions, giving the coincidence distribution as 
a function of — Xf. The data of both figures fit well to gaussian 
functions. We also see that the tails of the gaussians go to zero, 
indicating low background counts outside the coincidence region. 

Next we performed scanning measurements, this time using the 
SLM as described in the previous section. Figure 5 shows several two- 
dimensional scans, obtained by scanning in the vertical direction of 
the signal and idler photons. The scans were performed using phase 
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Figure 5 | Experimental results for 2D scans. The top row corresponds to raw data, giving distributions N. The bottom row corresponds to corrected data 
C. The splitting of the elliptical shapes seen in the smaller slit widths is probably due to the mapping of the intensity distribution of the pump-laser 
beam, which happens when the coincidences are measured in the near field. 
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Figure 6 | Experimental results for the marginal distributions of the signal photon. Red circles correspond to raw data Ng, and blue circles to data 
processed according to the procedure described in the paper, giving C^. 



slits of widths: 40, 36 (not shown), 30, 20, 15 and 10 pixels, as 
depicted in the different graphs. In all measurements (with SLM 
and with scanning detectors) accidental -coincidence count rates, 
defined as Q^c ^ c^C/At, were subtracted. In this formula, and C/ 
are the detector single-photon count rates and At is the temporal 
coincidence window. The upper row shows the raw uncorrected data, 
giving the raw distributions N. The usual elliptical pattern showing 
spatial correlations can be seen. As the slit width decreases in size, an 
increasingly prominent background contribution can be observed, 
which appears as the cross-shaped pattern covering the entire scan- 
ning region. It can be seen that the contrast decreases as the slit size 
decreases. The lower row of Figure 5 shows the data when the cor- 
rection procedure described above has been applied, giving the cor- 
rected distributions C. We can observe that the background counts 
are greatly reduced. The correction procedure is less effective for 
smaller slit sizes, where the signal-to-noise ratio is smaller. As the 
phase noise in the SLM is time dependent^" within the millisecond 
range, this noise is more effective when the photon flux is small. 

To better observe the reduction in background counts, in Figure 6 
we show one- dimensional marginal distributions for only the signal 
photon as a function of the slit position. Red circles correspond to 
uncorrected data and blue squares to corrected data C^. Gaussian 
curve fits serve as a guide to the eye. We observe that the inhomo- 
genous background contribution, represented by the flat tails of the 
gaussian distributions, is nearly eliminated in the corrected measure- 
ments. These plots can be compared to Fig. 4a), which was obtained 
by scanning the detector Di directly. Using the distribution shown in 
Fig. 4a), we calculate the variance of the distributions Cs(Xs) on the 
SLM, giving (Ax^) = (5274 + 170) pixels^. Converting the width of 
the mechanical slit in terms of the dimensions of the SLM pixels, and 
taking into account the magniflcation of the optical system, the result 
obtained is nearly 13 pixels. Therefore, we compare it to the variance 



of the corrected distribution obtained with a 15 pixel width phase slit, 
obtaining (5771 ± 380) pixels^, showing compatibility to within 
experimental error. Calculating the variance directly from the uncor- 
rected data gives a value that is incorrect by several orders of 
magnitude. 

To quantify the improvement in background counts, we define the 
contrast C = Cmax/(^taib where Cmax is the count rate at the maximum 
of the gaussian, and Cfau is the nearly constant count rate at the flat tail 

80 , , , , , , , , 



70 




Slit size (pixels) 

Figure 7 | Contrast as a function of the width of the phase slit, for 
uncorrected measurements Ns (red circles) and corrected measurements 
Cs (blue squares). 
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Figure 8 | Marginal distributions for the anti-diagonal coordinate Xg — Xj. Red circles correspond to raw data, giving distribution N_, and blue squares 
correspond to corrected data, giving distribution C_. 



of the gaussian distribution. In Figure 7 we plot the contrast for both 
uncorrected and corrected measurements as a function of the slit 
size, and do a linear fit to the data. The ratio of the slopes of the lines 
is about 15.75, showing that we obtained an improvement in the 
contrast of more than an order of magnitude for the corrected data. 

By summing the coincidence distributions (N or C) over the diag- 
onal variable Xs + Xf, we have the marginal coincidence distributions 
N_ and C_ for the anti-diagonal variable Xs — X/, as shown in Fig. 8. 
Red circles correspond to raw data N_ and blue squares to corrected 
data C_. We can see that the corrected data more closely approx- 
imates the distribution C{Xs — x^) shown in Fig. 4b). 

Our method is the simplest correction, in the sense that we assume 
that the value of ps {pd is the same for all pixels of the signal (idler) 
photon. We believe that more sophisticated methods are also viable, 
if one takes into account more technical details about the SLM. 

Conclusion 

In conclusion, we have presented a method for performing scanning 
measurements of spatial correlations of photons using an imperfect 
phase-only SLM. Exploiting the polarization dependence of the 
device, we use polarization interference measurements to turn a 
phase mask into an amplitude mask. We present a model that takes 
into account inhomogenous background counts introduced by the 
SLM and other imperfections in the system. We observe a consid- 
erable increase in contrast. We show that our method is useful for 
characterizing and studying spatial correlations of photon pairs. 
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